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Abstract 

We study the coherent transport modeled by continuous-time quantum walks, fo- 
cussing on hierarchical structures. For these we use Husimi cacti, lattices dual to 
the dendrimers. We find that the transport depends strongly on the initial site of 
the excitation. For systems of sizes N < 21, we find that processes which start at 
central sites are nearly recurrent. Furthermore, we compare the classical limiting 
probability distribution to the long time average of the quantum mechanical transi- 
tion probability which shows characteristic patterns. We succeed in finding a good 
lower bound for the (space) average of the quantum mechanical probability to be 
still or again at the initial site. 
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Introduction 

The dynamics of excitons is a problem of long standing in molecular and 
polymer physics (0). Thus, the incoherent exciton transport can be efficiently 
modeled by random walks, see, for instance, (JUIH); then, the transport follows 
a master (rate) equation and the underlying topology is fundamental for the 
dynamics (0). Studying the coherent transport, we model the process using 
Schrodinger's equation, which is mathematically closely related to the master 
equation, where the transfer rates enter through the connectivity matrix A of 
the structure (jH H; 0) . Moreover, also the elements of the secular matrix in 
Huckel's theory are given by A (J8|). 

Hyperbranched macromolecules have attracted a lot of attention in recent 
years. In this respect dendrimers have served as a prime example (0). Among 
a series of very interesting and crossdisciplinary applications like drug deliv- 
ery the role of dendrimers as light harvesting antennae has been investigated 
fjiol : 11; 0; EH). The details of the transport of an optical excitation over 



these structures depend on the localized states used. In one picture one may 
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envisage the excitation to occupy the branching points of the dendrimer, see, 
for instance, ESQ- However, as shown in fl5l). it may happen that the 
exciton occupies preferentially the bonds between the branching points. Then, 
the essential underlying structure is different and is given by sites localized at 
the mid-points of the bonds. We exemplify the situation in Fig. 1(a), starting 
from a dendrimer of generation 2 (open circles) and indicating the mid-points 
of the bonds by filled circles. Connecting neighboring filled circles by new 
bonds, one is led to the dual lattice of the dendrimer, which is called a Husimi 
cactus. Figure 1 shows three finite Husimi cacti of sizes N = 9, 21, and 45. 



(a) 





Fig. 1. Finite Husimi cacti (filled circles) of size (a) N = 9, (b) N = 21, and 
(c) N = 45. (a) also shows with dashed lines and open circles the corresponding 
dendrimer. 



Coherent dynamics modelled by quantum walks 

The quantum mechanical extension of a continuous-time random walk (CTRW) 
on a network (graph) of connected nodes is called a continuous-time quantum 
walk (CTQW). It is obtained by identifying the Hamiltonian of the system 
with the (classical) transfer matrix, H = — T, see e.g. (0; 0) (we will set H = 1 
and m = 1 in the following). The transfer matrix of the walk, T = (TJj), is 
related to the connectivity matrix A of the graph by T = — 7A, where for 
simplicity we assume the transmission rates 7 of all bonds to be equal and 
set 7 = 1 in the following. The matrix A has as non-diagonal elements 
the values —1 if nodes i and j of the graph are connected by a bond and 
otherwise. The diagonal elements An of A equal the number of bonds, fi, 
which exit from node i. 

The basis vectors \j) associated with the nodes j span the whole accessible 
Hilbert space to be considered here. The time evolution of a state \j) starting 
at time t is given by \j;t) = XJ(t,t )\j), where XJ(t,t ) = exp[— iH(t — t )] 
is the quantum mechanical time evolution operator. The transition amplitude 
Qtk,j{t) from state \j) at time to state \k) at time t reads then ockj(t) = 
(k\U(t,Q)\j) and obeys Schrodinger's equation. Denoting the orthonormalized 
eigenstates of the Hamiltonian H = — T by \q n ) (such that J2 n \ln){<in\ = 1) 
and the corresponding eigenvalues by A n , the quantum mechanical transition 
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probability (TP) is 



7r feJ (0 = \a kJ (t)\' 



-i\ n t 



Qn)(q n \j) 



(1) 



Note that for the corresponding classical probabilities J2kPkj{t) = 1 holds, 
with Pk,j(t) = Z)n(^l e_An< lQ , n) (Qn\j), whereas quantum mechanically we have 

E fc K,'(t)| 2 = i. 

Transition probabilities 

In the following we present the TPs 7r k ,j{t) obtained from diagonalizing the 
Hamiltonian H by using the standard software package MATLAB. For the 
finite Husimi cactus consisting of 21 nodes, as depicted in Fig. 1, we find 
numerically that the TPs are nearly periodic when the initial excitation is 
placed on one of the (symmetrically equivalent) nodes 1, 2, or 3 of the inner 
triangle. Thus at time t = 35.551846 one has 7Ti,i(35.551846) = 0.998761, 
when starting from 71"! i(0) = 1. In fact, the following expression holds exactly 
for all k sites of the cactus with N = 21: 



a f ] cos(Ajt) 



+ 



5> (fc) sin(A ; t) 



U=l 



(2) 



Here Ai = 0, A 2 = 3, A 3 = 3 + y/2, A 4 = 3 - ^2, A 5 = 3 + 2y/2, A 6 = 3 - 2y/2, 

(k) 

and the a\ follow from the corresponding eigenvectors. We further remark 
that the N = 21 cactus has 8 distinct eigenvalues, but that not all of them 
enter 7r M (t). 

(a) 1, — , — , — , — , — , — , — n— , (b) 
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Fig. 2. Transition probabilities TTjj(t) of the transport over finite Husimi cacti of 
sizes N = 21 (full lines) and N = 45 (dashed lines) for (a) j = 1 and (b) j = N. 

Figure 2 shows the probability to be at the initial node after time t for two 
different sizes iV and two different initial conditions. For TV = 21, when the 
excitation starts at the core (j = 1), the TP ni t i(t) is nearly recurrent. This is 
not anymore the case for iV = 45. Moreover, when the excitation starts at the 
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periphery (J = N), the TPs 7r/v,jv(i) are even more irregular. Now the iTk,N(t) 
for N = 21 with j = N indeed also contain the two eigenvalues A7 = 3 + y/E 
and As = 3 — \/6 in addition to those given after Eq. (2), a fact which increases 
the irregularity of the temporal behavior of the 7Tk,N(t)- 
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Fig. 3. Classical transition probabilities Pjj(t) of the transport over finite Husimi 
cacti of sizes TV = 21 (full lines) and N = 45 (dashed lines) for (a) j = 1 and (b) 
j = N. 

For comparison we also show in Fig. 3 the classical transition probabilities 
Pjj(t) for the same initial conditions as in the previous CTQW. Also here, 
there is a difference between the starting points. If the excitation starts at the 
central node 1, the equipartitioned limiting probability 1/N is reached much 
faster than when starting at the periphery. 

Long time limit 

The quantum mechanical time evolution is symmetric to inversion. This pre- 
vents 7ifc j(t) from having a definite limit for t — > 00. In order to compare the 
classical long time probability with the quantum mechanical one, one usually 
uses the limiting probability (LP) distribution (Jlfih 



1 r T 

Xk,j = lim - / dt TTkjit), 

T^oa I Jo 



(3) 



which can be rewritten by using the orthonormalized eigenstates of the Hamil- 
tonian, \q n ), as (0; E3) 



XK 



£5a„ 

n,m 



,x m (k\Qn)(q n \j)(j\qm)(q m \k). 



(4) 



Some eigenvalues of H might be degenerate, so that the sum in Eq. (4) can 
contain terms belonging to different eigenstates \q n ) and \q m ). 

In Fig. 4 we present the LPs for two sizes of Husimi cacti, iV = 21 and iV = 45. 
The LP distributions are self-similar generation after generation. Furthermore, 
for each size there are LPs having the same value, i.e. Xk,j — Xi,j- We collect 
LPs of the same value into clusters. Depending on where the excitation starts, 
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Fig. 4. Limiting probabilites for finite Husimi cacti of sizes N = 21 (left) and N = 4:5 
(right). 

the clustering is different. For instance, for N = 21, with an excitation starting 
at node 1, there are 6 different clusters, namely node 1 forms one cluster 
and the nodes {2, 3}, {4, 5}, {6, 7, 8, 9}, {10, 11, 12, 13}, and {14, • • ■ , 21} the 
remaining 5 clusters. When the excitation starts at another node, the clusters 
are formed by different nodes, as can be infered from Fig. 4. 

These results are closely related to our findings for the coherent transport 
over dendrimers (|l4h. For dendrimers the grouping of sites into clusters also 
depends on the starting site. Furthermore, there exists a lower bound for the 
LPs; one has namely Xk,j > 1/N 2 for all nodes k and j (|14h . This lower bound 
is related to the fact that in our cases A has exactly one vanishing eigenvalue, 
A = 0. 

Averaged probabilities 

As discussed above, starting at a central node of the N = 21 cactus, say node 
1, leads to nearly periodic 7Tfci(t). However, for larger cacti and/or different 
initial conditions this does not have to be the case. 

Classically one has a very simple expression for the probability to be still or 
again at the initially excited node, spatially averaged over all nodes. Then one 
finds 0) 

TV 1 N 



Ptf) = M^PiAt) 6XP ( ~ Xnt )- 

JV j=l JV n=l 



(5) 



This result is quite remarkable, since it depends only on the eigenvalue spec- 
trum of the connectivity matrix, but not on the eigenvectors. 

Quantum mechanically, we obtain a lower bound by using the Cauchy-Schwarz 
inequality, i.e., (14), 



1 N 1 r 

7F ( t ) = M ^ ^i-A 1 ) > eX P - *( A « _ X rn)t 
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In analogy to the classical case, the lower bound in Eq. (6) depends only on 
the eigenvalues and not on the eigenvectors. Note that for a CTQW on a 
simple regular network with periodic boundary conditions the lower bound 
is exact. We sketch the proof for a ID regular network of length N. There, 
the transition amplitudes to be still or again at node j after time t read 



a 



3,3 



(t) = N' 1 E„exp(-iA n £), see Eq. (16) of (|17|) with A n = 2-2 cos{2irn/N). 



The average 7r(t) is then given by 



W(t) 
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N ^ 
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n,m 



i{K ~ A m )t • (7) 



This result also holds for hypercubic lattices in higher, (^-dimensional spaces, 
when the problem separates in every direction; then one has ocjj(t) = [djj(t)] d 
(for the d — 2 case see (0)). 
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Fig. 5. Comparison of the averaged probability to be still or again at the initial node 
with the lower bound given in Eq. (6) for finite Husimi cacti of sizes (a) N = 21 
and (b) N = 45. 

In Fig. 5 we compare the exact value of w(t) with the lower bound given in 
Eq. (6). We use cacti of two different sizes, namely iV = 21 and iV = 45. 
Now, the lower bound fluctuates more than the exact curve, but it reproduces 
quite well the overall behavior of the extrema of the exact curve. Especially 
the strong maxima of the exact value are quantitatively reproduced. This is 
in accordance with previous studies on exciton transport over square lattices 
and over dendrimers (|14l). 



Conclusions 

In conclusion, we have modeled the coherent dynamics on finite Husimi cacti 
by continuous-time quantum walks. The transport is only determined by the 
topology, i.e. by the connectivity matrix, of the cacti. For the N = 21 cactus 
the dynamics is nearly periodic when one of the central nodes is initially 
excited. For larger structures and/or different initial conditions we observe 
only partial recurrences. To compare these results to those of the classical 
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(incoherent) case, we calculated the long time average of the TPs. Depending 
on the initial conditions, these show characteristic patterns, by which different 
nodes displaying the same limiting probabilities may be grouped into clusters. 
Furthermore, we calculated a lower bound for the average probability to be 
still or again at the initial node after some time t. This lower bound depends 
only on the eigenvalues of A and agrees quite well with the exact value. 
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